Numerical manifold method for dynamic nonlinear analysis of saturated porous media     

H. W. Zhang  L. Zhou 《国际地质力学数值与分析法杂志》2006,30(9):927-951

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions     

Kostantinos Papadakis  Constantine Goudas 《Astrophysics and Space Science》2006,305(2):99-124

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.  相似文献   

On the Assessment of Failure in Slope Stability Analysis by the Finite Element Method   总被引：3，自引：0，他引：3  

H. Zheng  D. F. Liu  C. G. Li 《Rock Mechanics and Rock Engineering》2008,41(4):629-639

  相似文献   

多传感器观测下带乘性噪声系统的逆向滤波与反褶积融合算法     

褚东升  李玉全  梁猛  张玲 《中国海洋大学学报(自然科学版)》2002,32(5):837-842

针对多传感器观测环境下带乘性噪声系统的逆向最优滤波与反褶积融合估计问题 ,本文提出了 1种基于极大似然准则的最优融合算法。该算法中各单传感器间并行计算 ,并且融合中心与单传感器处理中心间无反向通讯 ,因而执行效率较高。仿真表明 ,该融合算法产生的逆向滤波与反褶积比单传感器处理结果有较明显提高  相似文献   

高压水射流与潜水员水下作业安全研究     

张国光 《海洋工程》2006,24(3):113-118

从高压水射流技术在水下工程中应用及高压水射流的作用效应和伤害特征着手,分析影响水下高压水射流作业安全的因素,总结各国制订高压水射流作业的安全标准规程现状,提出制订我国水下高压水射流作业安全规程若干认识和设想。  相似文献   

波浪物理模型试验中风速比尺确定方法初探     

唐筱宁  夏运强  杨洪旗  贾友海 《海岸工程》2006,25(1):1-5

在波浪模型试验中研究风对构筑物的作用时,会涉及到风速比尺的确定问题。风对构筑物的作用力包括惯性力和粘滞力,这2种力的大小与建筑物的形态有关,对有的构筑物的作用力以惯性力为主,有的则以粘滞力为主。所以,在波浪物理模型试验中不能简单地采用重力相似准则也不能采用粘滞力相似准则来确定风速比尺。根据大量的试验研究结果,提出了采用不同模型比尺的船舶模型试验来确定风速比尺的方法,为波浪模型试验中风场模拟和数据处理提供了依据。  相似文献   

Change of multifractal thermal characteristics in the western Philippine sea upper layer during internal wave-soliton propagation     

Peter C. Chu  Chung-Ping Hsieh 《Journal of Oceanography》2007,63(6):927-939

The upper layer (above 140 m depth) temperature in the western Philippine Sea near Taiwan was sampled using a coastal monitoring buoy (CMB) with 15 attached thermistors during July 28–August 7, 2005. The data were collected every 10 min at 1, 3, 5, 10, 15, and 20 m using the CMB sensors, and every 15 sec at 15 different depths between 25 m and 140 m. Internal waves and solitons were identified from the time-depth plot of the temperature field. Without the internal waves and solitons, the power spectra, structure functions, and singular measures (representing the intermittency) of temperature field satisfy the power law with multi-scale characteristics at all depths. The internal waves do not change the basic characteristics of the multifractal structure. However, the internal solitons change the power exponent of the power spectra drastically, especially in the low wave number domain; they also break down the power law of the structure function and increase the intermittency parameter. The physical mechanisms causing these different effects need to be explored further.  相似文献   

南海北部冬季和夏季逆风流机制初探Ⅱ.季风逆风流产生机制     

袁叔尧  邓九仔 《热带海洋学报》1997,(3)

分析讨论了季风风应力、大陆坡地形及底摩擦在产生季风逆风流的必要条件和间接逆风流诊断判据中的作用；应用季风逆风流必要条件和间接逆风流诊断判据，解释冬季风和夏季风逆风流是如何产生的。结果表明：季风风应力是产生季风逆风流的主导因素；冬季风风应力、大陆坡地形及底摩擦三者联合作用导致表层海水在大陆坡上产生辐合生成高水位带和高动力高度带，在大陆被北侧产生NE向的冬季逆风流，南侧产生SW向顺风流；夏季风风应力、大陆坡地形及底摩擦三者联合作用导致表层海水在大陆坡上产生辐散生成低水位带和低动力高度带，在大陆被北侧产生SW向的夏季逆风流，南侧产生NE向的顺风流；冬季风盛行期间，风致经巴士海峡流入南海的黑潮水，将加速冬季逆风流的形成，加大冬季逆风流的强度；夏季风盛行期间，风应力的作用使巴上海峡以东的黑潮水不能进入南海，即使别的原因令巴上海峡以东的黑潮水流入南海，但高温、高盐的黑潮水对夏季逆风流具有阻扼作用。  相似文献   

利用Jason-1卫星确定全球平均海面高     

赵倩  黄舒青  朱明磊 《测绘信息与工程》2008,33(1):4-6

介绍了卫星测高数据处理的基本方法,分析研究了卫星测高的主要误差,提出了Jason-1的数据编辑准则,并用其测高数据确定了全球平均海面高,与CLS01模型进行了比较分析.  相似文献   

基于权变理论的滑坡灾害监测预报思路探讨     

胡华  俞登荣 《中国地质灾害与防治学报》2008,19(3):38-41

文章简要评述了地质灾害基础理论与应用技术发展现状、滑坡灾害多种监测预报判据的利弊。利用综合信息处理决策方法,提出了基于权变理论的滑坡灾害监测预报新思路。分析了滑坡成灾的权变特征、环境因素和决策因素,建立了滑坡灾害预报决策概念模型。进一步探讨了在预报决策中应遵循的动态性及满意性原则,为提高地质灾害监测预报理论的科学性提供了新的理论依据与技术途径。  相似文献